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文章 编号 :1000- 4939 (2023)06- 1428- 09 


未 知 测量 噪声 下 的 结构 多 类 型 啊 应 重 构 


史 鹏 程 , 彭 珍 瑞 


(兰州 交通 大 学 机 电工 程 学 院 ,730070 兰州 ) 


摘 要 :针对 使 用 传统 卡尔 曼 滤波 ( Kalman filter, KF) 算法 进行 结构 响应 重 构 时 测量 噪声 方差 估计 
及 多 类 型 响应 重 构 难 的 问题 ,提出 了 一 种 融合 噪声 方差 估计 和 多 类 型 响应 重 构 的 方法 。 首 先 在 移 
动 窗口 卡尔 受 滤波 算法 中 引入 近似 计算 的 方法 ,提出 递 推 移动 窗口 卡尔 曼 滤 波 算 法 (recursive mov- 
ing-window Kalman filter,RMWKF) 实时 北 推 估计 测量 噪声 方差 ,然后 基于 有 限 位 置 的 应 变 片 和 加 速 
度 传感器 的 测量 数据 ,并 结合 模 态 法 对 结构 各 个 位 置 的 应 变 、 位 移 、 速 度 以 及 加 速度 响应 进行 重 构 。 
最 后 对 二 维 析 架 和 号 居 梁 分 别 进 行 数值 模拟 和 试验 分 析 , 验 证 了 所 提 方 法 的 有 效 性 。 试 验 结果 表 
明 ,所 提 方 法 能 对 测量 噪声 方差 进行 有 效 地 估计 ,能 实现 有 效 的 多 类 型 结构 响应 重 构 ,同时 减 小 了 
数据 储存 空间 和 计算 量 。 
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Multi-type structural responses reconstruction 
with unknown measurement noise 


SHI Pengcheng, PENG Zhenrui 


(School of Mechanical Engineering ,Lanzhou Jiaotong University,730070 Lanzhou ,China ) 


Abstract :In order to solve the problem that it is difficult to estimate the variance of the measurement noise 
and reconstruct the multi-type response when the traditional Kalman filter ( KF) algorithm is used to 
reconstruct the structural response ,a method combining the noise variance estimation and the multi-type 
response reconstruction is proposed. First, the approximate calculation method is introduced into the 
moving window Kalman filter ( MWKF ) , and the recursive moving window Kalman filter( RMWKF) is 
proposed to estimate the real-time noise variance. Then ,the strain ,displacement ,velocity and acceleration 
responses of each position of the structure are reconstructed by using the measured data of strain gauges 
and acceleration sensors at limited positions combined with the modal method. Finally, the proposed 
method are verified by numerical simulation and experimental analysis of two-dimensional truss and 
cantilever beam respectively. The results show that the proposed method can estimate the variance of 


measurement noise and can reconstruct multi-type structural responses effectively. Moreover, the data 
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storage space and computation are reduced. 
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各 类 结构 动态 响应 的 测量 数据 是 进行 结构 健康 
监测 的 基础 。 目 前 各 种 机 械 、 土 木工 程 结构 都 日 益 
复杂 多 样 ,在 结构 健康 监测 中 所 需要 安装 传感器 的 
数量 随 之 增加 。 然 而 ,由 于 传感器 安装 位 置 限制 以 
及 经 济 等 原因 ,不 可 能 在 所 有 需要 位 置 都 安装 传 感 
器 ,或 由 于 部 分 传感器 长 期 服役 后 损坏 ,从 而 导致 进 
行 结构 健康 监测 所 需 的 数据 丢失 。 因 此 基于 有 限 的 
传感器 测量 数据 重 构 出 其 他 位 置 各 类 动态 响应 的 结 
构 响 应 重 构 技 术 就 具有 十 分 重要 的 研究 意义 1) 。 

目前 针对 结构 响应 重 构 技术 的 研究 可 依据 是 否 
基于 精确 的 有 限 元 模型 分 为 2 类 :确定 性 结构 响应 
重 构 和 不 确定 性 结构 响应 重 构 。 确 定 的 结构 响应 重 
构 方法 基于 精确 的 有 限 元 模型 ,建立 已 知 测量 响应 
与 待 重 构 响 应 之 间 的 关系 ,实现 对 未 知 响应 的 重 构 ， 
具体 方法 包括 基于 模 态 分 析 " 1 、 基 于 传递 率 矩 
阵 ' ”基于 经 验 模 态 分 解 ”""、 基 于 马尔 科 夫 
链 中 中 等 方法 。 但 对 于 一 些 较 为 复杂 的 结构 ,要 获 
得 其 精确 的 有 限 元 模型 是 较为 困难 的 , 且 其 结构 参 
数 在 长 时 间 的 服役 后 必然 会 有 所 改变 。 若 有 限 元 模 
型 的 精度 得 不 到 保证 , 则 必然 影响 结构 响应 重 构 的 
精确 性 。 在 不 确定 性 结构 响应 重 构 中 ,不 必 基 于 精 
确 的 有 限 元 模型 ,同时 考虑 到 数据 采集 过 程 中 存在 
噪声 的 影响 ,可 以 采用 卡尔 曼 滤波 .扩展 卡尔 曼 滤波 
等 卡尔 曼 滤波 类 算法 进行 结构 响应 重 构 。 卡 尔 曼 滤 
波 (Kalman filter,KF) 算 法 是 一 种 采用 递 推 形式 的 时 
域 滤波 方法 ,能够 有 效 处 理 模型 误差 和 测量 误差 , 实 
现 对 系统 状态 的 最 优 估计 。 张 笑 华 等 "基于 KF 算 
法 对 二 维 检 架 的 应 变 和 位 移 响 应 进行 了 重 构 , 同 时 
对 位 移 和 应 变 两 种 传感器 的 数量 和 位 置 进行 了 优 
化 。XU 等 外 使 用 KF 对 大 跨度 悬索桥 进行 应 变 、 位 
移 .加 速度 的 响应 重 构 ,并 结合 逐步 消去 法 对 多 类 型 
传感器 进行 优化 布置 。 任 鹏 等 使 用 分 别 增 广 卡 
尔 曼 滤波 方法 和 传递 率 方法 在 未 知 外 部 激励 的 情况 
下 对 一 钢 检 架 进 行 了 应 变 响 应 重 构 , 并 对 比 了 2 种 
方法 的 优 劣 。 董 康 立 等 591 基于 类 卡尔 曼 滤 波 算法 
验证 了 对 外 部 激励 进行 重 构 的 有 效 性 ,并 以 响应 重 
构 精 度 为 目标 使 用 董 火 虫 算法 和 逐步 消去 法 进行 传 
感 器 优化 布置 。 

以 上 不 确定 性 的 响应 重 构 中 ,大 多 基于 卡尔 曼 


滤波 类 算法 ,在 卡尔 曼 滤 波 类 算法 中 ,系统 噪声 方差 
抢 阵 测量 噪声 方差 矩阵 以 及 系统 的 初始 状态 值 等 
参数 对 于 滤波 结果 非常 重要 ,使 用 不 准确 的 噪声 方 
差 矩 阵 会 降低 滤波 的 精度 ,甚至 可 能 导致 滤波 结 
严重 发 散 … 。 在 实际 测量 过 程 中 噪声 方差 矩阵 不 
是 一 个 固定 值 , 而 是 随时 间 不 断 变化 的 ,但 现 有 
的 基于 卡尔 曼 滤波 类 算法 的 结构 响应 重 构 中 大 多 将 
上 述 测量 噪声 方差 矩阵 设置 为 一 个 定 值 ,这 显然 是 
不 准确 的 。 为 此 ,LAI 等 ”采用 了 移动 窗口 扩展 卡 
尔 曼 滤波 算法 对 噪声 方差 矩阵 进行 在 线 估计 ,用 于 
一 个 6 层 楼 结构 的 损伤 识别 ,验证 了 此 方法 对 噪声 
的 强 适应 性 和 和 鲁 棒 性 。ZHANG 等 2 在 此 基础 上 使 
用 了 移动 窗口 卡尔 曼 滤 波 算法 (moving-window Kal- 
man filter, MWKF) 对 一 简 支 梁 的 应 变 和 位 移 信 号 进 
行 重 构 。 

本 研究 在 上 述 研究 的 基础 上 ,针对 在 使 用 KF 
算法 进行 响应 重 构 时 , 因 测量 噪声 方差 无 法 准确 获 
得 ,从 而 造成 重 构 精 度 降低 的 问题 ,引入 递 推 移动 窗 
口 卡尔 曼 滤 波 法 (recursive moving-window Kalman 
filter, RMWKF) ,实现 了 对 测量 噪声 方差 的 在 线 递 推 
估计 ,与 MWKF 相 比 减 小 了 数据 的 储存 空间 和 计算 
量 。 并 结合 模 态 法 ,在 仅 使 用 工程 中 常用 的 应 变 片 
和 加 速度 传感器 的 情况 下 ,基于 有 限 点 的 测量 数据 
对 结构 各 处 的 应 变 、 位 移 .速度 和 加 速度 响应 进行 重 
构 ,最 后 对 二 维 棉 架 和 悬臂 深 分 别 进行 数值 模拟 和 
试验 验证 。 


1 理论 基础 


1.1 动力 学 方程 离散 化 


结构 的 二 阶 运动 方程 用 模 态 坐标 表示 为 
dg +2€ wd + wg = BD'Lu (1) 
其 中 :9 为 模 态 坐标 ;é 为 阻尼 矩阵; oo 为 模 态 频率 
和 矩阵; 为 外 部 激励 向 量 ; 民 为 激励 映射 矩阵 ;更 为 
模 态 位 移 振 型 矩阵 。 
将 式 (1) 写成 状态 空间 方程 并 离散 化 


| =Ax,+Bu,+w: 


(2) 


z= Cx,+Du,+yv 
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其 中 : x; 、zs 和 到 分 别 为 离散 后 的 状态 向 量 、 观 测 
向 量 和 外 部 激励 向 量 ; 4 和 B 分别 为 离散 化 后 的 状 
态 和 矩阵 和 输入 和 矩阵; C 和 D 分 别 为 输出 矩阵 和 直接 
传输 矩阵 ; w 和 六 分别 为 由 于 模型 的 不 确定 性 造成 
的 系统 噪声 和 测量 误差 造成 的 测量 噪声 ,分别 假设 
为 均值 为 零 的 高 斯 白 噪 声 ; 方 差 矩 阵 分 别 为 2 和 R。 
式 (2) 中 的 各 项 具体 表示 为 


人 
Pe 中 | 
4 = as 本 (4) 
0 I 加 0 
| | UD[ | 2 
三 0 
C = 0 oD (6) 
-Do -2Ge om 
0 
D | 0 (7) 
DDL 


其 中 ,五 =T7G,7 为 应 变 -位 移 转 换 和 矩阵 。 
1.2 应 变 -位 移 转换 矩阵 


结构 的 位 移 可 以 通过 其 模 态 坐标 和 位 移 振 型 的 
全 加 表示 为 


d= 之 9(z)9， (8) 
其 中 :mm 为 位 移 模 态 阶 数 ;g, 为 第 > 阶 模 态 坐标 ; 
2 (Xx) 为 第 r 阶 模 态 位 移 振 型 。 
位 移 和 应 变 之 间 的 关系 可 以 表示 为 


5 = (9) 
~ Ox 
其 中 ,s, 为 x 方向 的 应 变 , 将 式 (8) 代 入 (9) 得 
$= (Dp)g,) (10) 
结构 应 变 可 以 表示 为 
$= 3 pg, (11) 
W(x) = rz) (12) 


中 ,更 (x) 为 应 变 模 态 。 将 式 (12) 代 入 式 (11) 中 
得 


5 = 工序 (04 (13) 
由 式 (8) 可 得 所 有 已 知 节点 位 移 可 以 表示 为 


和 力学 学 报 


第 40 卷 
d, = 94 (14) 
其 中 : d, 为 所 有 已 知 的 节点 位 移 ; 2. 为 已 知 测 点 的 
横 态 振 型 矩阵 ,由 于 大 部 分 情况 下 gp。 都 不 是 方 阵 ， 
要 求解 4 不 能 直接 求 逆 , 因此 需 两 边 同 时 左 乘 


T 4 日 


9 和合 


ps d,=9p, pd (15) 
式 (15) 中 gp,"y, 为 方 阵 ,两 边 同 时 左 乘 (oro,) 


得 


4=(2o。2,) pd, (16) 
将 式 (16) 代 入 式 (11) 得 
_9 T 二 
S$, 二 rn Pe 9,) 9。 ad, (17) 
_9 T A 
T= Pp 90,) Ou (18) 


其 中 ,gp, 为 全 部 节点 的 模 态 位 移 振 型 矩阵 。 由 此 可 
得 到 应 变 位 移 转 换 和 矩阵 7。 


1.3 ”KF 算法 


卡尔 曼 滤 波 算 法 是 一 种 针对 线性 系统 ,基于 状 
态 空间 方程 ,通过 输入 观测 数据 ,对 系统 状态 进行 最 
优 估计 的 算法 。 其 大 致 可 以 分 为 预测 和 更 新 2 个 步 
又 ,预测 部 分 生成 移 验 状态 估计 和 先 验 估计 误差 协 
方差 ;更 新 部 分 结合 移 验 估计 和 卡尔 曼 增益 ,生成 后 
验 状 态 估计 和 后 验 估计 误差 协 方差 。 


预测 : 
xis = Axntt+Bu, (19) 
Pi,=AP,,A' +0 (20) 

更 新 : 
Ki = PC ICP CT+RI! (21) 

ee 三 Rind + 用 [Zi 一 Ge — Du,,!] 

22) 
Pi = [TK CIP (23) 


其 中 : x 为 先 验 状态 估计 ;Pi ,11 为 先 验 估计 误差 
方差 窍 阵 ;KK , ;为 卡尔 曼 增 益 矩 阵 ; sR 为 后 验 
状态 估计 ;Pi ,14; 1 为 后 验 估计 误差 方差 矩阵 ;R 和 
2 分 别 为 测量 噪声 方差 矩阵 和 系统 噪声 方差 矩阵 。 


1.4 MWKF 算法 


由 于 在 传统 卡尔 曼 滤 波 算 法 中 测量 噪声 方差 矩 
阵 灵 常常 被 定义 为 一 个 常量 ,但 实际 过 程 中 测量 噪 
声 方 差 是 时 变 的 ,使 用 不 准确 的 测量 噪声 方差 会 导 
致 滤波 精度 的 降低 甚至 滤波 的 严重 发 散 。 为 了 计算 
出 实际 测量 过 程 中 的 测量 噪声 方差 ,文献 [20] 中 首 


il 
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次 引入 MWKF 算法 对 测量 噪声 协 方 差 进行 估计 。 

由 于 在 响应 重 构 中 真实 的 测量 响应 是 难以 获得 
的 ,从 而 各 个 时 刻 真 实 的 测量 噪声 &; 也 无 法 准确 佑 
计 。 因 此 采用 测量 值 的 加 权 平 均 来 近似 获得 真实 
响应 吓 ] 。 

ZU = OQ Zh + Oy Zh 十 Qs ZX (24) 

其 中 : a as as 为 加 权 系 数 ; zi 表示 真实 测量 响应 ， 
从 而 测量 噪声 应 为 


EE 二- 放 访 (25) 


R, 分 别 可 以 写 为 ' 


E(z,) > (26) 
R, 二 i 1 lS -天 (si)) (2 -ECs))] 
(27) 


其 中 ,WN 为 统计 样本 个 数 。 
在 计算 最 终 的 测量 噪声 方差 时 将 多 个 移动 窗口 
进行 加 权 平 均 处 理 , 即 
Ri; = yi Ra + yy Ri + + yi Rs (28) 
其 中 : y; 为 加 权 系 数 ; i 为 采样 窗口 个 数 ,计算 如 下 
1/tr(R,) 
/rR ) + /tr(Rs) + + 1/tr(R,) 
(29) 


Yi 


其 中 , tr(: ) 为 计算 矩阵 的 迹 。 
1.5 RMWKF 算法 


在 上 述 MWKF 算法 中 ,由 于 求解 最 终 的 测量 吧 
声 方差 时 ,采用 了 加 权 平均 处 理 ,从 而 导致 数据 的 储 
存 空间 和 计算 量 大 大 增加 。 为 了 减 小 数据 的 存储 和 
计算 量 ,引入 了 递 推 估计 的 方法 ,实现 了 对 于 测量 吕 
声 方差 的 递 推 估计 。 将 式 (27) 拆 分 为 
Es-Ee))(s -Es))"]+ 


1 
(es) (esl :0) 


根据 式 (27) 可 知 Ri 为 


天 -1 


R = HTS- EB(e)) (as ~- BCs))"] 
(31) 

式 (31) 可 以 近似 写 为 

Re = 3 Bs Ele)) (6 ~ ECs))"] 
(32) 
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综合 式 (30) 和 式 (32) , Ri 的 递 推 佑 计 为 
N-2 
R, = Nw- Re + 


a 


由 以 上 对 测量 噪声 方差 矩阵 Ri 的 递 推 过 程 可 
以 看 出 ,RMWKF 相 比 于 MWKF'” 实现 了 对 测量 噪 
声 方差 的 在 线 递 推 估计 ,kk 时 刻 的 测量 噪声 方差 R, 
可 由 -1 时 刻 的 测量 噪声 方差 R, 1 和 时 刻 的 噪 
声 值 及 均值 递 推 得 到 。 在 保持 与 MWKF 算法 重 构 
精度 相近 的 前 提 下 , RMWKF 只 需 计 算 时 刻 测量 
噪声 及 其 均值 ,结合 -1 时 刻 的 测量 噪声 方差 即 可 
算得 大 时 刻 的 测量 噪声 方差 ,但 MWKF 则 需要 对 窗 
口内 w 个 数据 进行 计算 处 理 。 显 然 , 相 比 于 
MWKF,RMWKF 具有 计算 量 小 和 占用 储存 空间 小 的 
优势 。 


1.6 结构 响应 重 构 


运用 上 述 递 推 移动 窗口 卡尔 曼 滤波 算法 进行 状 
态 的 预测 和 更 新 ,可 以 得 到 响应 重 构 方程 , 即 
| = Cx, +Du, 


(34) 
zi = C, Xx, +D, uu 
其 中 : z; 和 必 分 别 为 需要 重 构 位 置 的 真实 计算 响应 
和 重 构 响应 ,其 中 包括 应 变 响应 速度 响应 和 加 速度 
响应 ; C; 和 D; 分 别 表示 需要 重 构 位 置 的 输出 矩阵 
和 传输 矩阵 ,分 别 为 C 和 忆 的 子 矩 阵 ,下 标 * 为 选择 
的 模 态 ,在 实际 使 用 中 只 需 考虑 对 响应 贡献 较 大 的 
前 几 阶 模 态 。 求 解 出 结构 的 应 变 响应 后 ,可 通过 位 
移 与 应 变 之 间 的 关系 ,对 结构 的 位 移 进 行 重 构 。 根 
据 模 态 站 加 原理 ,结构 位 移 和 应 变 分 别 可 以 表示 为 
d,s = 更 wm dn (35) 
,a (36) 
其 中 :n, 表 示 结 构 自 由 度数 ;n, 表 示 结 构 单元 数 ;m 
表示 所 选取 的 模 态 阶 数 ;结合 式 (35) 和 (36) 可 以 直 
接 通 过 应 变 响应 求 出 结构 的 位 移 响应 为 
di = ,a Wn (37) 
其 中 ,上 标 + 表示 求 伪 道 。 
由 于 结构 应 变 和 加 速度 两 种 响应 在 数量 级 上 相 
差 很 大 ,这 会 引起 输出 矩阵 C 的 高 度 病态 ,导致 卡尔 
曼 滤波 结果 的 不 准确 ,因此 要 对 输出 矩阵 进行 规则 
化 处 理 后 代入 卡尔 曼 滤波 算法 进行 计算 
C= (R) CC" (38) 
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2 数值 模拟 


采用 二 维 杭 架 数值 模拟 整个 响应 重 构 过 程 , 以 
验证 所 提 方 法 的 有 效 性 。 如 图 1 所 示 ,该 枯 架 结构 
有 14 节点 ,25 个 单元 和 25 个 自由 度 。 弹 性 模 量 、 密 
度 和 泊 松 比分 别 为 210 GPa、7 850 kg/m 和 0.3。 选 
取 前 六 阶 模 态 振 型 为 目标 模 态 进行 响应 重 构 ,前 六 
阶 振 型 频率 分 别 为 17. 72 .58. 89 .85. 48 .113. 67、 
169. 98 .220.22 Hz。 其 中 前 六 阶 振 型 的 具体 形式 分 
别 为 第 一 、 第 四 和 第 六 阶 为 对 称 垂 向 弯曲 ;第 二 和 第 


由 


图 2 给 出 了 第 11 单元 的 应 变 以 及 4 节点 的 了 
方向 的 加 速度 响应 使 用 RMWKF 算法 计算 的 测量 品 
声 方差 的 变化 情况 ,从 图 中 可 以 看 出 传感器 的 测量 


第 40 卷 


五 阶 为 反对 称 垂 向 弯曲 ;第 三 阶 为 横向 拉 伸 。 在 第 
4 市 点 的 了 方向 施加 随机 激励 ,激励 的 频带 宽度 为 
10 ~250 Hz, 最 大 值 为 35 N。 按 照 文献 [16] 中 所 提 
出 的 逐步 消去 法 优化 布置 传感器 的 布置 ,在 第 10、 
11 .15 和 第 16 单元 布置 4 个 应 变 片 ,第 3 和 第 4 个 
节点 的 了 方向 布置 2 个 加 速度 传感器 。 利 用 状态 空 
间 方 程式 (34) 算 出 布置 传感器 处 的 真实 计算 响应 ， 
并 在 计算 的 应 变 和 加 速度 值 中 分 别 加 入 均值 为 0， 


方差 为 1 x10 一 和 2.56 的 高 斯 白 噪 声 作为 测量 
响应 。 


4 
上 m 


1.5m 
图 1 二 维 检 架 结构 


Fig.1 Two-dimensional truss structure 


计算 ,实验 测试 环境 为 :AMD Ryzen 53500U 处 理 右 ， 
8G 内 存 , 基 于 MATLAB 2018a 环境 采用 M 语言 编写 
算法 ,窗口 长 度 V 设 定 为 20, 运 行 时 间 结 果 如 表 1 


噪声 方差 随 着 时 间 不 断 变化 ,计算 的 噪声 方差 在 整 
个 过 程 中 围绕 预 设 的 真实 值 上 下 波动 ,由 此 可 知 使 
用 RMWKF 算法 估计 的 测量 噪声 方差 是 接近 预 设 值 
的 ,可 以 用 于 后 续 的 响应 重 构 计 算 。 


0 D2 SROLLE2 LI M6 S20 
ts 


(a) 11 单 元 应 变 

7 

6 

所 

| 

& 

2F， 

1 

0 02 "04 0:06 0:30 10 2 4 J 2.0 


(b) 1 加 速度 
图 2 传感器 测量 噪声 方差 
Fig.2 Covariance of estimated measurement noise 
为 了 验证 RMWKF 算法 在 实际 计算 测量 噪声 方 
差 时 , 相 比 于 MWKF 算法 是 否 拥有 计算 量 小 和 占用 
储存 空间 小 的 优势 。 在 保证 计算 精度 的 前 提 下 ,分 
别 对 于 不 同 测量 长 度 的 测量 数据 进行 测量 噪声 方差 


所 示 。 可 以 看 出 RMWKF 算法 对 于 同一 长 度 的 测量 
数据 ,计算 时 间 仅 为 MWKF 算法 的 三 分 之 一 左右 ， 
进一步 验证 了 RMWKF 具有 计算 量 小 和 占用 储存 空 
间 小 的 优势 。 

表 1 2 种 算法 运行 时 间 对 比 (单位 :s) 


Tab.1 Comparison of running time of two algorithms (unit:s) 
要 数据 长 度 
算法 
4 000 8 000 20 000 40 000 50 000 


RMWKF 0. 109 0.213 0.460 0.767 0.874 


MWKF 0.304 0.553 1.216 2.041 2.553 


将 RMWKF 算法 的 重 构 效果 与 相同 条 件 下 使 用 
KF 和 MWKF 算法 进行 重 构 的 效果 进行 对 比 ,使 用 
相对 百分比 误差 Pu 来 衡量 各 自重 构 的 精度 , 且 
std(z —2z) 
i x 100% (39) 

图 3 是 分 别 采用 KF.MWKF 以 及 RMWKF 进行 
响应 重 构 的 各 个 单元 以 及 节点 的 各 类 响应 的 Pi 的 
对 比 结果 ,在 测量 噪声 方差 未 知 的 情况 下 采用 RM- 
WKF 响应 重 构 的 Pi 要略 低 于 采用 KF 响应 重 构 的 
Pi ,与 采用 MWKF 进行 重 构 的 相对 百分比 误差 相 
近 。 图 4 所 示 是 结构 各 处 的 响应 ,采用 RMWKF 进行 
响应 重 构 的 各 类 响应 时 程 与 计算 响应 时 程 吻合 良好 ， 
可 以 用 于 测量 噪声 方差 未 知情 况 下 的 响应 重 构 计 算 。 


Par = 
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(b) 5 节点 7 向 位 移 
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(d) 9 节点 7 向 加 速度 


图 4 重 构 响 应 与 计算 响应 对 比 


Fig.4 Comparison of reconstruction and computational responses 


3 试验 验证 


使 用 图 5 所 示 甚 臂 梁 作为 试验 对 象 , 进一步 对 
以 上 方法 进行 试验 验证 。 悬 臂 梁 尺寸 为 600 mm x 
50 mm x7. 8 mm, 弹性 模 量 为 128 GPa, 密度 为 
6 440 kgjm , 泊 松 比 为 0.3。 将 悬臂 粱 分 为 12 个 单 
元 ,13 个 节点 和 24 个 自由 度 。 选 取 前 四 阶 模 态 振 
型 为 目标 模 态 进行 响应 重 构 ,前 四 阶 振 型 频率 分 别 
为 1$.73 .96.38 .274.92 .3$43. 18 Hz。 其 中 前 四 阶 振 
型 的 具体 形式 分 别 为 :第 一 阶 为 垂 向 单 向 弯曲 ,第 二 
阶 为 垂 向 双向 弯曲 ,第 三 阶 为 反对 称 垂 向 弯曲 ,第 四 


阶 为 对 称 垂 向 弯曲 。 在 节点 9 的 了 方向 施加 激励 ， 
激励 时 程 曲线 如 图 6 所 示 。 按 照 与 数值 模拟 中 相同 
的 逐步 消去 法 确定 传感器 的 位 置 , 如 图 7 所 示 ,在 第 
5 个 单元 布置 一 个 应 变 片 ,第 9、 第 10 以 及 第 13 个 
节点 布置 3 个 加 速度 传感器 。 并 在 已 布置 传感器 的 
基础 上 ,在 第 2 个 单元 和 第 7 个 节点 分 别 布置 一 个 
应 变 片 和 一 个 加 速度 传感器 作对 比 以 检验 重 构 效 
果 。 试验 所 采用 的 应 变 片 型 号 为 BFH120-3AA; 应 
变调 理 器 型 号 为 INV1861A; 加 速度 传感器 型 号 为 
ICP 型 的 INV9821 ;网 络 分 布 式 采集 仪 为 北京 东方 振 
动 和 噪声 技术 研究 所 的 INV3062-C2。 
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时 测量 噪声 方差 的 变化 曲线 。 
到 5 甚 臂 梁 实物 
Fig.5 Cantilever beam 0 Qn O10 015 020 0.30 0.35 0.40 0.45 0.50 
ts 
, (a) 5 单元 应 变 

之 -2 
沁 -4 
渗 -6 
8 
一 10 


0 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 
ts 


多 6 激励 


Fig.6 Excitation 


图 应 变 传感器 
全 加 速度 传感器 


图 7 传感器 布置 位 置 图 
Fig.7 Location of sensor arrangement 

为 保证 动 应 变 的 测量 精度 ,应 注意 以 下 方面 : 
将 胶粘剂 均匀 地 涂抹 在 完好 的 应 变 片 上 ,将 应 变 
片 粘贴 在 经 过 单 向 打磨 光滑 并 清洗 的 被 测 结构 表 
面 ,并 对 其 表面 进行 涂 覆 保护 ; @) 采用 温度 补偿 应 
变 片 来 降低 温度 效应 对 测试 结果 的 影响 ;@ 根 据 被 
测 结构 的 受 力 确 定 粘贴 位 置 和 电 桥 的 组 成 方式 , 消 
除 附 加 载荷 , 测 取 单一 载荷 ;由 确保 信号 采集 盒 接 
地 ,应 变调 理 絮 调 平衡 ,测试 系统 调 零 , 再 进行 动 应 
变 测 量 。 

虽然 传感器 附加 质量 与 悬臂 梁 整体 质量 相 比 较 
小 ,但 是 在 试验 过 程 中 仍然 会 对 梁 的 局 部 模 态 造成 
一 定 的 影响 。 以 下 两 个 方面 尽量 减 小 传感器 附加 质 
量 带 来 的 影响 。 人 测量 时 尽量 选择 质量 更 轻 的 传 感 
器 ,并 用 动力 学 修改 预测 软件 来 消除 附加 质量 带 来 
的 影响 , 即 在 传感器 测量 位 置 添加 一 个 相同 质量 的 
“ 负 质 量 " 来 消除 此 影响 ;@@ 由 于 采用 基于 卡尔 曼 滤 
波 的 不 确定 性 响应 重 构 方法 ,可 允许 有 一 定 的 建 模 
误差 ,附加 质量 不 会 对 最 终 重 构 结 果 造 成 较 大 的 
影响 。 

在 根据 式 (24) 计算 真实 响应 时 ,加 权 系 数 a 、 
0 分 别 取 1/6 .2/3 .176; 在 计算 测量 噪声 方差 的 
数据 统计 样本 个 数 NN 取 20。 图 8 给 出 了 位 于 5 单 
元 和 9 节点 的 一 个 应 变 片 和 一 个 加 速度 传感器 测量 


l/s 
(b) 9 节点 7 向 加 速度 


图 8 传 感 费 测量 噪声 方差 

Fig.8 Covariance of estimated measurement noise 

基于 RMWKF 对 测量 噪声 方差 估计 后 进行 响应 
重 构 计 算 , 并 分 别 与 使 用 KF 和 MWKF 重 构 的 结 
对 比 。 表 2 列 出 了 部 分 单元 和 节点 重 构 信 号 的 相对 
百分比 误差 。 


表 2 相对 百分比 误差 


Tab.2 Relative percentage errors 


Prg/% 
单元 /节点 号 
KF MWEKF RMWKF 

应 变 

2 0.52 0.23 0.22 

9 3.33 1.25 1.24 
位 移 

7 0.30 0.11 0.15 

11 0.22 0.07 0.07 
速度 

7 1.43 0.61 0.61 

11 0.71 0.30 0.30 
加 速度 

7 5.82 2.09 2.09 

11 3.81 1.38 1.38 


从 表 2 中 可 以 看 出 采用 RMWKF 进行 响应 重 构 
的 结果 的 精度 与 使 用 MWKF 的 精度 类 似 , 同 时 优 于 
KF 的 精度 ,证 实 了 采用 移动 窗口 法 对 于 测量 噪声 方 
差 进 行 实时 估计 的 重要 性 。 图 9 和 图 10 分 别 给 出 
了 重 构 响 应 与 测量 响应 和 计算 响应 的 对 比 ,使 用 
RMWKF 结合 模 态 法 的 响应 重 构 方 法 可 以 仅 在 使 用 
应 变 片 和 加 速度 传感器 的 情况 下 ,实现 对 于 结构 各 
处 的 应 变 \ 位 移 、 速 度 和 加 速度 信息 的 重 构 。 重 构 的 


投稿 网 站 :http://cjam. xjtu. edu. en 微 信 公 众 号 :应 用 力学 学 报 


第 6 期 


言 号 数据 与 传感器 测量 数据 真实 计算 数据 吻合 良 
好 ,可 以 重 构 出 未 布置 传 感 带 处 的 结构 的 各 类 真实 
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(b) 7 节点 7 向 加 速度 


9 重 构 响 应 与 测量 响应 对 比 


Fig.9 Comparison of the reconstruction and measured responses 
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4 结 论 


1) 在 仅 使 用 应 变 片 和 加 速度 传感器 的 情况 下 ， 


使 用 RMWKF 结合 模 态 法 可 以 重 构 出 结构 各 点 的 应 
变 .位移 .速度 以 及 加 速度 响应 , 且 各 类 重 构 的 响应 
与 计算 以 及 测量 响应 均 吻 合 良好 。 

2) 基 于 RMWKF 算法 可 对 测量 噪声 方差 进行 实 
时 估计 ,估计 结果 基本 接近 预 设 值 ,与 固定 的 测量 噪 
声 方 差 相 比 ,使 用 时 变 的 测量 噪声 方差 更 符合 工程 
实际 , 且 相 比 于 MWKF 算法 ,其 具有 计算 量 小 和 占 
用 储存 空间 小 的 优势 。 

3) 相 比 于 使 用 KF 进行 响应 重 构 , RMWKF 可 在 
测量 噪声 方差 未 知 的 情况 下 进行 结构 响应 重 构 , 且 
数值 算 例 和 试验 皆 验 证 了 使 用 RMWKF 的 重 构 精 度 
与 使 用 MWKF 重 构 的 精度 相近 , 相 比 使 用 KF 重 构 
的 精度 要 更 高 。 

本 研究 的 数值 算 例 和 实验 对 象 结构 相对 简单 ， 
作者 正在 此 基础 上 进一步 拓展 研究 ,将 本 研究 所 述 
的 方法 应 用 于 大 型 复杂 结构 中 ,以 验证 应 用 于 工程 
实际 中 的 可 靠 性 。 
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